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ABSTRACT 



We present Rossi X-ray Timing Explorer observations of the X-ray pulsar 
SMC X-1. The source is highly variable on short time scales (< 1 h), exhibiting 
apparent X-ray flares occupying a signiflcant fraction (~3 %) of the total observ- 
ing time, with a recurrence time of ~100 s. The flares seem to occur over all 
binary orbital phases, and correlate with the overall variability in the light curve. 
We find a total of 323 discrete flares which have a mean full width half maximum 
of ~18 s. The detailed properties of SMC X-1 do not vary significantly between 
the flares and the normal state, suggesting that the flare may be an extension of 
the normal state persistent emission with increased accretion rates. The flares 
resemble Type II X-ray bursts from GRO J1744-28. We discuss the origin of the 
SMC X-1 flares in terms of a viscous instability near the inner edge of the accre- 
tion disk around a weakly magnetized X-ray pulsar, and flnd this is consistent 
with the interpretation that SMC X-1 is in fact an intermediate-stage source like 
GRO J1744-28. 

Subject headings: accretion, accretion disks — pulsars: individual (SMC X-1) — 
stars: neutron — X-rays: bursts — X-rays: stars 



1. Introduction 

Neutron star X-ray binaries are generally categorized into two groups: low mass X-ray 
binaries (LMXBs) and X-ray pulsars. The surface magnetic fleld of the central neutron star 
in an LMXB is thought to be ~10^ G. The mass accretion with this magnetic fleld is most 
likely spherical, so that no signiflcant inhomogeneity in the X-ray emission over the neutron 
star surface is expected - i.e., no persistent coherent pulsations are observed. The strong 
magnetic fleld, ~10^^ G, of an X-ray pulsar, on the other hand, can funnel the accretion 
matter onto the magnetic pole, which makes the central neutron star appear as a pulsar. 
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Of particular interest are the so-called "intermediate-stage sources" speculated to lie 
between the LMXBs and X-ray pulsars, including "the Rapid Burster" (MXB 1730-355), 
GRO J1744-28 ("the Bursting Pulsar"), and SAX J1808.4-3658 ("the accreting millisecond 
pulsar") (Lewin et al. 1976; Fishman et al. 1995; in 't Zand et al. 1998). Both the 
Rapid Burster and GRO J1744-28 exhibit Type 11 X-ray bursts; however, only the former 
shows Type 1 X-ray bursts while only the latter has apparent coherent pulsations (Lewin et 
al. 1996). SAX J1808. 4-3658, on the other hand, shows both Type I bursts and coherent 
pulsations (Wijnands & van der Klis 1998; Chakarbarty & Morgan 1998), but not Type II 
bursts. The magnetic field strengths of these sources have been inferred to be ~10^~^^ G, 
between those of LMXBs and X-ray pulsars. 

Another possible intermediate-stage source is the X-ray pulsar SMC X-1, which has 
similar properties to GRO J1744-28, including its fast spin period (~0.72 s for SMC X-1; 
~0.47 s for GRO J1744-28), steady spin-up, and inferred magnetic field (~10^^ G ) (Bildsten 
& Brown 1997; Li & van den Heuvel 1997). In addition, once SMC X-1 was observed with 
an X-ray burst that resembles Type 11 bursts (Angelini, Stella, & White 1991). It may be 
possible, therefore, that SMC X-1 and GRO J1744-28 belong to a distinctive group of X-ray 
binaries, "bursting pulsars", which show both coherent pulsations and Type II X-ray bursts 
(Li & van den Heuvel 1997). To investigate this important possibility, we analyze all publicly 
available RXTE data for SMC X-1, searching for phenomena that may be related to X-ray 
bursts. We report that SMC X-1 in fact exhibits active flares resembhng Type II bursts from 
GRO J1744-28. 



2. Data Analysis and a Flcire Search 

We analyzed all publicly available RXTE Proportional Counter Array observations to- 
ward SMC X-1. Photon arrival times from the Good Xenon mode were transformed to the 
solar system barycenter using the JPL DE400 ephemeris. The Very Large Event Models 
were used to subtract backgrounds, and only the Standard 2 data obtained from the top 
xenon layers of Proportional Counter Units 0, 1, and 4 were considered for spectral analysis. 
Both the data within 30 minutes after passages through South Atlantic Anomaly and/or 
with a high (> 0.1) electron ratio were ignored. 

The data toward SMC X-1 are occasionally contaminated by outbursts of the nearby 
(~27' away) transient source XTE JOlll.2-7317 (Chakrabarty et al. 1998). The contam- 
ination was easily identified by the source's ~31-s pulsations, with powers reach up to > 
5000 in Leahy-normalized PDSs. We searched for weak contamination in PDSs as narrow 
statistically significant (> 10 in Leahy-nomalized PDSs) peaks at ~0.032 Hz (and/or at its 
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harmonic frequencies). We excluded any data contaminated by XTE JOlll. 2-7317 from our 

analyses, as well as data containing dip- like features (e.g., Her X-1; Moon & Eikenberry 
2001) or close to the eclipse (i.e., binary orbital phase between 0.9 and 0.1) of SMC X-1. 
Wc obtained a total of 150 data segments with an average length of ~1360 s, giving a net 
observational duration of ~204 ks. For spectral analysis, we considered only the spectrum 
between 2.5-25 keV range due to a PCA responsivity problem (R. Remillard 2001, private 
communication), and assumed a systematic uncertainty of 1 %. 

In our analysis, we define "fiarcs" to be the part of a light curve that has three or more 
consecutive data bins with photon counts larger than a threshold value of 3a Poisson noise 
above the mean photon count in a given light curve. We analyzed all the 150 light curves 
as follows searching for the flares. First, we binned each fight curve to 4-s time resolution 
and made a flare fist. We extended the search with 2- and 8-s time resolutions, excluding 
the flares already found with different time resolutions. We found a total of 323 flares, and 
fit them a Gaussian function to estimate width and peak intensity. We obtained ~18 s and 
~15 s for the mean and rms standard deviation of FWHM; the reduced (= xl) of fits 
range from 0.5 to 1.9. The total duration above half-maximum is ~5.8 ks, indicating that 
SMC X-1 spends ~3 % of its time on flares. 



3. Flare Examples and Correlations with Other Parameters 

The rms variability of the 150 light curves ranges between 6 and 16 %, with a mean of 
~11 %. Figure 1 presents three light curves (and their PDSs) with very low (7 %), average (11 
%), and very high (16 %) rms variability as examples representing three different variability 
levels. The mean photon count rates are 120.5 ± 8.5, 116.1 ± 12.9, and 146.4 ± 23.3, for 
Figure la, lb, and Ic, respectively. No ffare is found in Figure la, which has a low rms 
variability; three flares are found Figure lb, which has an average rms variability. Strong 
flaring activity is well illustrated in Figure Ic, which has a high rms variability, with several 
apparent flares lasting for a few tens of seconds. At the flare peak, the photon count rate 
rises up to ~2.5 times of that outside the flares. The inset shows a Gaussian fit to the 
intense flare at t ~ 900 s, with an FWHM estimated to be ~23 s. Leahy-normalized PDSs 
of the three light curves (Figure Id, Ic, and If) show strong QPO-like peaks aroimd 10 
mHz, as well as the peaks caused by the source's coherent pulsations at z/ ~ 1.4 Hz (and 
their harmonics), indicating the existence of ~100-s aperiodic variability independent of the 
flaring activity. The fractional rms (FRMS) amplitudes in the 2-50 mHz range are 7 ± 1, 9 
± 1, and 14 ± 1 % for Id, le, and If, respectively. 

Figure 2 presents the pulse profiles of the three light curves, all showing the double- 
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peaked, smooth profile typical of SMC X-1 (e.g., Levine et al. 1993). They have the same 
pulsed fractions of ~39 %, with their second peaks at phase 0.56 with respect to the first 
ones at phase 0. The ratio of the second peak to the first peak is ~0.78 for Figure 2a, while 
it is ~0.96 for 2b and 2c. No significant variation in the pulse peak ratio has been found 
within a given light curve, regardless whether it is obtained inside or outside flares. 

Figure 3 compares the phase-averaged softness ratio, defined to be the ratio of the soft 
X-ray (2-5 keV) photon count rates to those of the hard X-ray (5-13 keV), of the three light 
curves in Figure 1 with the total photon count rates at 2-25 keV with 4-s resolution. The 
average softness ratio of the three light curves is invariant: 0.47 ± 0.03, 0.47 ± 0.02 , and 
0.46 ± 0.04 for Figure la, lb, and Ic, respectively. This contrasts with the softness ratio of 
the LMC X-4 flares (Moon, Eikenberry, & Wasserman 2002), which shows a strong hnear 
correlation with the total photon count rate (inset in Figure 3). 

We examined the spectral invariance of SMC X-1 over the flaring activity indicated by 
the constant softness ratio distribution more thoroughly via fitting the 32-s spectra obtained 
from the largest peak in each of the three hght curves (i.e., the peak at f ~ 900 s in Figure 
la, at i ~ 1050 s in lb, and at t ~ 1110 s in Ic) to a model spectrum. The model spectrum 
consists of a power- law component with a high-energy cutoff (for non-thermal magnetospheric 
emission) and a Gaussian component (for iron line emission), together with a component 
for photoelectric absorption by intervening interstellar matter. We fixed the central energy 
of the Gaussian component to be 6.7 keV based on previous results (Angelini et al. 1991). 
Figure 4 compares the observed spectra with the best-fit model spectra, and Table 1 lists the 
best-fit parameters and xt of fits. Although all parameters are poorly constrained, the 
power-law index does not change significantly over the flaring activity - consistent with the 
softness ratio distribution (Figure 3). The difference in TVh may be caused by the motion of 
the precessing, tilted accretion disk suggested to be responsible for the super-orbital motion 
of SMC X-1, but the low energy hmit of the spectral flts (i.e., 2.5 keV) makes it difficult 
to constrain A^h propertly, because the photoelectric absorption by intervening interstellar 
matter is expected to be most signiflcant in the soft energy band. 

In order to perform statistical analyses, we calculated the flare fraction, which we deflne 
to be the ratio of the integrated time that SMC X-1 is flaring (i.e., within FWHMs of the 
Gaussian flts) to the total observing time of a given parameter, and investigate its correlation 
with the parameter. Figure 5 shows the distribution of the flare fraction as functions of the 
binary orbital phase (5a), rms variability of the light curve (5b), and the pulse peak ratio 
(5c). Some important results are worth noticing: the flare fraction (1) is larger than 1.4 % 
over all the orbital phases (with its minimum at the orbital phase of ~0.35), (2) shows a 
signiflcant variation from phase to phase, and (3) increases with the rms variability of the 
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light curve, as well as with the pulse peak ratio. The correlation between flare fraction and 
rms variability identified in Figure 5b remains very similar when we use the flare-subtracted 
rms variability. 

4. Discussion 

SMC X-1 shows active flares that occupy ~3 % of the total observing time. The flares 
seem to occur over all binary orbital phases, and the flaring activity is proportional to the 

rms variability of the light curve (Figure 5). Except for the small change in the pulse peak 
ratio (Figure 2 and 5), no significant change is found along with the flaring activity. All these 
suggest that the SMC X-1 flares may be simple extensions of the persistent emission of a 
normal state with increased accretion rates, but without significant changes in the geometry 
of the accretion flows and the magnetosphere. Under the hypothesis that the double peaks 
in the pulse profile (Figure 2) are due to the two magnetic poles of SMC X-1, one simple 
explanation for the change in the pulse peak ratio may be that the increase in the accretion 
rate onto the fainter pole is higher than that onto the brighter pole during the flares. The 
invariant pulse peak ratio in a given light curve (regardless of flares), together with the 
correlation between the flare-subtracted rms variability and the flare fraction, on the other 
hand, may indicate the existence of variability whose time scale is longer than the flare 
recurrence time, ~100 s. 

The SMC X-1 flares differ from Type I X-ray bursts from LMXBs for various reasons, 
including the shape of profiles and spectral properties. While the profile of Type I bursts 
shows an abrupt increase with an exponential decay in most cases, the SMC X-1 flares have 
symmetric Gaussian shapes. The X-ray spectrum of the SMC X-1 flare is far from thermal 
spectrum of Type I bursts (although the thermal spectrum of an X-ray pulsar is not well 
constrained, so it is not completely excluded that the SMC X-1 flare spectrum is thermal). 
In addition, the SMC X-1 spectrum does not show any apparent variation within a flare, 
while Type I burst shows spectral coohng as the burst continues. The changes in the softness 
ratio and pulse profile of the SMC X-1 flares differ from those of the LMC X-4 flares (Figure 
3; Moon et al. 2002), which is one of the most well known and regular flaring X-ray binaries. 

The SMC X-1 flares recall the Type II X-ray bursts found in GRO J1744-28 mainly 
due to the spectral invarince over the flaring activity. In fact, the burst spectrum of GRO 
J1744-28 was found to be very similar to that of SMC X-1, with a similar photon index of 
~1.2 and e-folding energy of ~14 keV (Lewin et al. 1996). One difference is that the SMC 
X-1 flares lack the post-flare dip which often follows the Type 11 bursts from GRO J 1744-28. 
The rather gradual rise of the SMC X-1 flares, however, may be responsible for it via offering 
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sufRcient time to replenish the material in the accretion disk. This is consistent with the 
interpretation that the SMC X-1 flares are simple extensions of a normal state with increased 
accretion rates, which is probably the most compelling argument for attributing the origin 
of the SMC X-1 flares to an accretion disk instability. Given the difference between GRO 
J 1744-28 and SMC X-1 (i.e., a transient low- mass system versus a persistent high- mass one), 
we consider that the magnetic field strength, suggested to be comparable for the two sources, 
is critical to understanding the Type II bursts from GRO J1744-28 and the X-ray flares from 
SMC X-1. 

Wc note that SMC X-1 may be capable of experiencing a viscous instability, namely the 
Lightman-Eardlcy instability (Lightman & Eardlcy 1974), due to its ~10^^ G cm^ magnetic 
moment (Li & van den Heuvel 1997). The radiation pressure in this case is comparable to 
the gas pressure around the inner disk radius (i.e., the magnetospheric radius), resulting 
in a viscous instability with shghtly increased mass accretion (e.g., Cannizzo 1996, 1997). 
Because the instability develops near the inner edge of the accretion disk, it has an advantage 
in explaining bursts/flares with a short recurrence time. The viscosity parameter (a) of the 
classical a-disk (Shakura & Sunyaev 1973) at the transition radius between the "inner region" 
and the "middle region" around an 1.4 Mq neutron star is a ~ 216 t~^^'^Mn, where t^is is the 
viscous time scale in seconds and Mn is the mass accretion rate in units of 10^^ g s^^. For a 
typical value M17 ~ 20 for SMC X-1 (e.g., Wojdowski et al. 2000), the viscosity parameter 
is q; ~ 0.14 and 4.4 when t^is is 1000 and 100 s, respectively. The value 4.4 for tyis = 100 s is 
somewhat larger than generally expected, a < 1. However, even larger values {a ~ 10-100) 
were obtained for the dwarf nova HT Gas, possibly due to the patch nature of the accretion 
disk (Vrielmann, Hessman, & Horne 2002). Alternately, the strong magnetic held of SMC 
X-1 may help increase the viscosity parameter (R. Lovelace 2002, private communication). 
We need detailed numerical studies to investigate this scenario more thoroughly. 



5. Summary and Conclusions 

Through the analysis of the all pubhcly available RXTE data toward the X-ray pulsar 
SMC X-1, we find that the source is highly variable on short time scales (< 1 h), and that 

Gaussian flares occur over all orbital phases with a ~100 s recurrence time scale. The flares 
occupy ~3 % of the total observing time, and the flaring activity is proportional to the 
overall variability of the source. While the pulse peak ratio shows a small systematic change 
along with the flaring activity, the PDSs, pulse proflles, softness ratios, and X-ray spectra 
during the flares are very similar to those outside the flares, indicating that the flares are 
probably extensions of a normal state just with increased accretion rates. This supports 
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the interpretation that the SMC X-1 flares have their origin in an accretion disk instabihty 
and the suggestion that it may belong to a distinctive group of "bursting pulsars" with the 
"bursting pulsar" GRO J1744-28, owing to its ~10^^ G surface magnetic field. A viscous 
instability near the inncr-cdgc of the accretion disk might be responsible for the SMC X-1 
flares, although detailed studies on this scenario need to be done in the future. 
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Table 1. Best-fit Parameters for Spectra Shown in Figure 4 



Parameter 


(a) 


(b) 


(c) 


A^H*(1022 cm-2) 


3.2(0.9) 


1.6(1.5) 


2.2(0.9) 


a* 


1.6(0.3) 


1.3(0.5) 


1.6(0.3) 


(keV)* 


17.2(7.1) 


6.9(3.2) 


14.3(5.2) 


(keV)* 


7.8(6.9) 


12.6(10.2) 


14.0(10.6) 




0.81 


1.3 


0.81 


Flux (10"^ ergs cm~^ s~^) 


1.7 


2.0 


2.8 



*7Vh and a are the hydrogen nuclei column density of the inter- 
vening matter and the index of the power law component. Ec and 
Ei represent the cutoff and e-folding energy of the high-energy 
cutoff component. 

**Each spectrum has 48 degrees of freedom. 



Note. — Energy range is 2.5-25 keV, and the 90 % uncertainty 
levels are quoted in the parentheses. 
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Fig. 1. — (a-c): Light curves (in the 2-25 keV range) with 4-s time resolution of three data 
segments that have ~7 % (a), ~11 % (b), and ~16 % (c) rms variabihties. The numbers 
in the upper-right corners represent the MJD of the beginning of each hght curve and its 
binary orbital phase. The inset in (a) compares the intense flare (crosses) at t ~ 900 s with 
a Gaussian fit (solid histogram). The arrows indicate the central positions of nine fiares. 
(d-f): Leahy- normalized PDSs of the three light curves in the left-hand side. 
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Fig. 2. — Pulse profiles of the three fight curves in Figure 2 - solid fine for 2a, dotted fine 
for 2b, and dashed line for 2c. The number of bin is 50. 
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Fig. 3. — Softness ratio of the three hght curves in Figure 1: triangle for la, crosses for lb, 
and diamonds for Ic. The softness ratio is a ratio of the soft X-ray (2-5 keV) photon count 
rates to those of hard X-ray (5-13 keV). The inset shows the same softness ratio of the LMC 
X-4 flares (Moon et al. 2002). All the softness ratios are obtained with 4-s resolution. 
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Fig. 4. — Comparison of the observed spectra (crosses) of the brightest peaks of the three 
Ught curves in Figure 1 (the peak at t ~ 900 s in la, t ~ 1050 s in lb, and t ~ 1110 s in Ic) 
with the spectra (sohd histogram) of the best-fit parameters. The bottom panels represent 
the residuals of the fittings. 
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Fig. 5. — The distributions of the flare fraction (sohd histograms) versus the binary orbital 
phase (a), the rms variabihty of the hght curve (b), and the pulse peak ratio (c). The error 
bars correspond to 1-a confidence level. The dashed histograms represent the integrated 
observing time (the right-hand ordinates) in the corresponding binary orbital phase (a), the 
rms variability (b), and the pulse peak ratio (c). 



